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The features of the emissions of light particles {LP), charged particles (CP), intermediate mass 
fragments (IMF) and the largest fragment (MAX) are investigated for 129 Xe as functions of 
temperature and "freeze-out" density in the frameworks of isospin dependent lattice gas model and 
classical molecular dynamics model. Definite turning points for the slopes of average multiplicity 
of LP, CP and IMF, and of the mean mass of the largest fragment (A max ) are shown around a 
liquid gas phase transition temperature and while the largest variances of the distributions of LP, 
CP, IMF, and MAX appear there. It indicates that the cluster emission rate can be taken as 
a probe of nuclear liquid gas phase transition and while the largest fluctuation is simultaneously 
accompanied at the point of the phase transition within lattice gas model and molecular dynamics 
model as the traditional thermodynamical theory requires when a phase transition occurs. 
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I. INTRODUCTION 

Phase transition and critical phenomenon is an exten- 
sive debating subject in nature sciences. Recently, the 
same concept was introduced into the microscopic sys- 
tems, such as in atomic cluster Q and nuclei The 
break-up of nuclei due to violent collisions into several 
intermediate mass fragments, can be viewed as critical 
phenomenon as observed in fluid, atomic, and other sys- 
tems. It prompts the possible signature on the liquid gas 
phase transition in the nuclear system. One fact is that 
the sudden opening of the multifragmentation Q| and 
vaporization Q channels can be interpreted as the signa- 
ture of the boundaries of phase mixture B. In addition, 
the fact that the nuclear caloric curve [[| 0, ||, |j| [H], 
in a certain excitation energy range shows a plateau gives 
a possible indication of a first order phase transition || 
as predicted in the framework of statistical equilibrium 
models [0|. On the other hand, the extraction of crit- 
ical exponents in the charge or mass distribution of the 
multifragmentation system |]l3| can be explained as an 
evidence of the phase transition. 

Several theoretical models have been developed to 
treat such a phase transition in nuclear disassembly. 
These models are percolation model jl4], [ll|, lattice gas 
model (LGM) and statistical multifragmentation model 
fl6f etc. In this work, what we are interested in is the 
lattice gas model. It is a simple short range interac- 
tion model |Tijl , but it can be successfully applied to 
nuclear systems with isospin symmetry and asymmetry 



|23|, |24|, [2§ |2§ |27|, |28|. LGM is car- 
ried assuming the system is in a "freeze-out" density pf 
with thermal equilibrium at temperature T. Previous 
calculations fl23l with LGM showed that there exists a 
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phase transition for the finite nuclear systems by study- 
ing the effective power law parameter (r) of cluster mass 
or charge distribution p9|, pCfl , their second moments (S2) 
|14j |, and the specific heat. More recently, we proposed 
two novel criteria, namely multiplicity information en- 
tropy (H) and nuclear Zipf's law to diagnose the onset 
of liquid gas phase transition in the framework of the 
isospin dependent LGM (I-LGM) and isospin dependent 
classical molecular dynamics (I-CMD) model J31| . 

In this work, we will show that the emission rate of 
clusters is a useful tool to diagnose the nuclear liquid gas 
phase transition, and while the largest fluctuation of clus- 
ter multiplicities is simultaneously revealed at the point 
of the phase transition by investigating the features of 
light particles, charged particles, intermediate mass frag- 
ments and the largest fragment of disassembling source 
in the frameworks of I-LGM and I-CMD. 

The paper is organized as following. Sec. II gives the 
descriptions of I-LGM and I-CMD. Results and discus- 
sions are shown in Sec. Ill where the multiplicities of 
cluster emissions, their slopes and their fluctuations are 
investigated. The influence of the "freeze-out" density 
on cluster emission is also presented in the framework of 
I-LGM, and comparisons of the results between I-LGM 
and I-CMD at a given "freeze-out" density is followed. 
Finally the conclusion is given in Sec. IV. 



II. DESCRIPTION OF MODELS 

A. Isospin dependent lattice gas model 

The lattice gas model of Lee and Yang |l7j , where the 
grand canonical partition function of a gas with one type 
of atoms is mapped into the canonical ensemble of an 
Ising model for spin 1/2 particles, has successfully de- 
scribed the liquid-gas phase transition for atomic system. 
The same model has already been applied to microscopic 
nuclear system. Detailed descriptions for the features of 
LGM can be found in some literatures, eg. papers of Pan 
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and Das Gupta et al. Jig, |l9|, of Campi and Krivine p0| , 
of Carmona, Richert and Tarancon |2l[| and of Chomaz 
and Gulminelli [^7) etc. Here we will adopt the model 
developed by Pan and Das Gupta (l8[ In order to 
better understand the context of this work, the models 
are described below. 

In the LGM, A (=N + Z) nucleons with an occupation 
number s (i.e. "spin") which is defined as s = 1 (-1) for 
a proton (neutron) or s — for a vacancy, are placed 
in the L sites of 3 dimension lattice. Nucleons in the 
nearest neighboring sites interact with an energy e SiS .. 
The hamiltonian of the system is written as 
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The interaction constant e SiS . is related to the binding 
energy of the nuclei. In order to incorporate the isospin 
effect in the lattice gas model, the short-range interac- 
tion constant e SiSj is chosen to be different between the 
nearest neighboring like nucleons and unlike nucleons: 



e pp = O.MeV, 
-5.33MeV, 



(2) 



which indicates the repulsion between the nearest- 
neighboring like nucleons and the attraction between 
the nearest- neighboring unlike nucleons. This kind of 
isospin dependent interaction incorporates, to a certain 
extent, Pauli exclusion principle and effectively avoids 
producing unreasonable clusters, such as di-proton and 
di-neutron clusters, etc. The disassembly of the system 
is to be calculated at an assumed "freeze-out" density 
Pf = (A/L)p , where po is the normal nucleon density, 
beyond pf nucleons are too far apart to interact. 

In the model, N + Z nucleons are put in L sites by 
Monte Carlo sampling using the canonical Metropolis al- 
gorithm |2(], [52). As pointed out in Ref. |21, 33 , however, 
one has to be careful treating the process of Metropolis 
sampling in order to satisfy the detailed balance prin- 
ciple and therefore warrant the correct equilibrium dis- 
tribution in the final state. Speaking in detail, in this 
work, first an initial configuration with N + Z nucleons 
is established. Second, for each event, sufficient num- 
ber of "spin" -exchange steps is tested, eg. 20000 steps 
to let the system generate states with a probability pro- 
portional to the Boltzmann probability distribution with 
Metropolis algorithm. In each "spin" -exchange step, a 
random trial change on the basis of the previous con- 
figuration is made. For instance, we choose a nucleon at 
random and attempt to exchange it with one of its neigh- 
boring nucleons or vacancies at random regardless of the 
sign of its " spin" (Kawasaki- like spin-exchange dynamics 
|54f ), then compute the change AE in the energy of the 
system due to the trial change. If AE is less than or 
equal to zero, accept the new configuration and repeat 
the next "spin" -exchange step. If AE is positive, com- 
pute the "transition probability" W = e~ AE / T and com- 
pare it with a random number r in the interval [0, 1]. If r 



< W, accept the new configuration; otherwise retain the 
previous configuration. 20000 "spin" -exchange steps are 
performed to assure we get the equilibrium state. Third, 
once the nucleons have been placed stably on the cubic 
lattice after 20000 " spin" -exchange steps for each event, 
their momenta are generated by Monte Carlo sampling 
of Maxwell Boltzmann distribution. Thus various ob- 
servables based on phase space can be calculated in a 
straightforward fashion for each event. One important 
point of such Monte-Carlo Metropolis computations is 
that the above "spin" -exchange approach between the 
nearest neighbors, independently of their "spin", is evi- 
denced to be satisfied by the detailed balance condition 
as noted in Ref. |2l], [53|. In other words, this sampling 
method will guarantee that the generated microscopic 
states form an equilibrium canonical ensemble. 

Once this is done the LGM immediately gives the clus- 
ter distribution by using the rule that two nucleons are 
part of the same cluster if their relative kinetic energy is 
insufficient to overcome the attractive bond : 



< 0. 



(3) 



This method has been proved to be similar to the 
so-called Coniglio-Klein's prescription in the condensed 
matter physics |35| and was shown to be valid in LGM 
g§ @ g(| ^, |f 



Isospin dependent classical molecular dynamics 
model 



Since the lattice gas model is a model of the near- 
est neighboring interaction, a long range Coulomb force 
is not amenable to lattice gas type calculation. Pan 
and Das Gupta provided a prescription, based on simple 
physical reasoning, to decide if two nucleons occupying 
neig hboring sites, form part of the same cluster or not 
fll9fl . They first try to map the lattice gas model cal- 
culation to a molecular dynamics type prediction, both 
first done without any Coulomb interaction. If the cal- 
culations match quiet faithfully then they can study the 
effects of Coulomb interaction by adding that to the 
molecular dynamics calculation. Here we adopt the same 
prescription to use the molecular dynamics and there- 
fore investigate Coulomb effect. The results and con- 
clusions can now be compared and checked between the 
LGM and the CMD. Obviously, here we do not do any 
ab initio molecular dynamics calculation but only use it 
for simulating the nuclear disassembly starting from a 
thermally equilibrated source which has been produced 
by the above I-LGM: i.e. the nucleons have been initial- 
ized at their lattice sites with the Metropolis sampling 
and have their initial momenta with Maxewell Boltzmann 
sampling. From this starting point we switch the calcu- 
lation to a classical molecular dynamics evolution under 
the influence of a chosen force. Noting that in this case 
Pf is, strictly speaking, not a "freeze-out" density for 
molecular dynamics calculation but merely defines the 
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starting point for time evolution. However since classical 
evolution of a many particle system is entirely determin- 
istic, the initialization does have in it all the information 
of the asympototic cluster distribution, we will continue 
to call pf the "freeze-out" density. 

I 



where rg = 1.842/m is the distance between the centers 
of two adjacent cubes. The parameters of the potentials 
axe p = 2, q = 1, a = 1.3, B = 0.924, and C = 1966 
MeV. With these parameters the potential is minimum 
at ro with the value -5.33 MeV, is zero when the nucleons 
are more than 1.3ro apart and becomes strongly repulsive 
when r is significantly less than ro. We now turn to the 
nuclear potential between like nucleons. Although we 
take e pp = e nn = in I-LGM, the fact that we don't put 
two like nucleons in the same cube would suggest that 
there is short range repulsion between them. We have 
taken the nuclear force between two like nucleons to be 
the same expressions as above plus 5.33 MeV up to r — 
1.842 fm and zero afterwards: 

v pp( r )i r < r °) = v n P {r) - v np (r ), 
v pp (r)(r>r0) = 0. (5) 

The system evolves with the above potential. The evolu- 
tion equations for each nucleon are, as usual, given by 

dpl/dt = -Y, j ^iS/ i v{r i j), 

dfl/dt = pi/m. (6) 

Numerically, the particles are propagated in the phase 
space by integrating Newton's equations of motion 
through a "leap-frog" algorithm. At asymptotic times, 
for instance, the original blob of matter has expanded 
to 64 times its volume in the initialization, the clusters 
are easily recognized: nucleons which stay together after 
arbitrarily long time are part of the same cluster. The 
observables based on cluster distribution in both mod- 
els can now be compared. While they can be also com- 
pared by switching on/off Coulomb interaction within the 
molecular dynamics. 

III. RESULTS AND DISCUSSIONS 

We choose a medium-size nucleus 129 Xe as an example. 
The input parameters are temperature T and "freeze- 
out" density pf in the model calculations. In this work 
T is mostly limited to the range of 3 to 7 MeV and the 
"freeze-out" density pf is chosen to be the range of 0.18po 



The form of the force in the CMD is chosen to be also 
isospin dependent in order to compare with the results 
of I-LGM. The potential for unlike nucleons is expressed 
as 



(4) 



to 0.60po- In most cases, pf is chosen to be 0.38 po, since 
the experimental data can be best fitted by a pf between 
0.3po and 0Ap o in previous LGM calculations (ll| [56), 
which corresponds to 7 x 7 x 7 cubic lattice for Xe. There 
is also good support from experiment that the value of 
Pf is significantly below 0.5po JUl- ln addition, 0.18po, 
corresponding to 9 x 9 x 9 cubic lattice and 0.60poj corre- 
sponding to 6 x 6 x 6 cubic lattice of " freeze-out" densities 
for 129 Xe are also taken to compare and check the results 
with different "freeze-out" density. 1000 events are sim- 
ulated for each combination of T and pf which ensures 
good statistics for results. 

A. I-LGM in different "freeze-out" density 

Fig.l shows that the mean multiplicities of emitted 
neutrons, protons, charged particles and intermediate 
mass fragments and the mean mass for the largest frag- 
ment evolve with temperature at different "freeze-out" 
densities in the I-LGM calculation. At a fixed "freeze- 
out" density, average neutron multiplicity (N n ), proton 
multiplicity (N p ), charged particle multiplicity (N cp ) and 
the largest fragment mass (A max ) display monotonously 
increasing or decreasing with temperature as expectedly. 
But the multiplicity (Nimf) of IMF shows a rise and 
fall with temperature |38| ^0| , when the system prob- 
ably crosses the phase transition boundary. With the 
decreasing of the "freeze-out" density, N p , N n , N cp and 
A m ax decreases since larger space separation among nu- 
cleons at smaller "freeze-out" density makes the clusters 
less bound and therefore decreases the sizes of free clus- 
ters and then increases the cluster multiplicities. The 
situation of Ni m f is slightly complicate, i.e. it increases 
with the decreasing of " freeze-out" density in lower tem- 
perature branch but it is contrary in high temperature 
branch. 



It seems to be impossible to discover the possibility of 
phase transition of nuclei if we only see these mean quan- 
tities as shown above (Ni m f is an exception). However, 
when we focus on their slopes to temperature: figure 2, 



v np (r)(r/r0<a) = C [B(r /r)r - (r /rf] exp( 1 — ), 
v np (r)(r/r > a) = 0. 

I 
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FIG. 1: Average multiplicity of the emitted neutrons N n (a), pro- 
tons N p (b), charged particles N cp (c) and intermediate mass frag- 
ment Ni m f (d), average mass of the largest fragment A max (e) as 
a function of temperature in different "freeze-out" density in the 
framework of I-LGM. 



the sharp changes are observed at the nearly same tem- 
perature at each fixed " freeze-out" density: for instance, 
4 MeV at 0.18p , 5 MeV at 0.38p and 6 MeV at 0.60p - 
At such a transition point, (1) the multiplicities of emit- 
ted clusters increase rapidly and after that the emission 
rate slows down; (2) the decrease of the largest fragment 
size reaches to a valley for such a finite system. Physi- 
cally the largest fragment is simply related to the order 
parameter pi - p g (the difference of density in nuclear 
"liquid" and "gas" phases). In the infinite matter, the 
infinite cluster exists only on the "liquid" side of critical 
point. In a finite matter, the largest cluster is present on 
both sides of phase transition point. In this calculation, 
a valley for the slope of A max to temperature may corre- 
spond to a sudden disappearance of infinite cluster ("bulk 
liquid") near the phase transition temperature. It is not 
occasional producing such waves of the slopes, it should 
reflect the onset of phase transition there. This idea can 
be supported by surveying the other phase transition ob- 
servables: such as the effective power law parameter, the 
Campi's second moment from the mass or charge dis- 
tribution of fragment and the information entropy H of 
event multiplicity distribution pq, pl[. Figure 3 depicts 



FIG. 2: Same as figure 1, but for their slopes with temperature. 

these results. The minimum of r and the maxima of S2 
and H appear around respective phase transition tem- 
peratures at different values of "freeze-out" density, i.e. 
about 4 - 4.25 MeV at 0.18p , 5 - 5.5 MeV at 0.38p and 
6 - 6.5 MeV at 0.60po- These temperatures are consistent 
with the ones extracted from the choppy position of the 
above slopes. It indicates that the above slopes (emis- 
sion rate) can be taken as a probe of a liquid gas phase 
transition of nuclei. 



Furthermore, the largest fluctuations of cluster mul- 
tiplicities are found around the phase transition point 
in the same calculation. Figure 4 illustrates that RMS 
width (a) of the multiplicity distributions of neutrons, 
protons, charged particles and intermediate mass frag- 
ments, and the distribution of the largest fragment 
masses. These variances generally show peaks at the 
same phase transition temperatures as the ones extracted 
from the above observables for each fixed "freeze-out" 
density. Noting that the fluctuation of A max is related to 
the compressibility of the system. These features are also 
consistent to the one of phase transition behaviors, i.e. 
there exists the largest fluctuation at the phase transition 
point |4l]| . This fluctuation represents internal feature of 
the disassembling system not numerical fluctuation. 
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FIG. 3: Phase transition observables: effective power law parame- 
ter t (a), the second moment 52 (b) of fragment mass distribution 
and the information entropy H (c) as a function of temperature in 
different "freeze-out" density in the framework of I-LGM. 



FIG. 4: Temperature dependent RMS widths (<r) of the distribu- 
tions of N n (a), N p (b), N cp (c), N im f (d) and Amax (e) in different 
"freeze-out" density in the framework of I-LGM. The symbols are 
the same as figure 1. 



We would like to point that the "freeze-out" density 
dependent phase transition temperature, extracting from 
t and 52, was also observed in a previous work |lS| ], but 
this pf dependent phenomenon will vanish when excita- 
tion energy is used as a variable Ejj. In other words, 
the excitation energy is perhaps a good correspondence 
of the critical temperature because of only one critical 
point and hence only one critical temperature for a sys- 
tem. 



B. Comparison of I-LGM to I-CMD in the same 
"freeze-out" density 

Considering the absence of long range Coulomb force 
in the LGM, we will adopt the CMD to investigate the 
Coulomb effect and check the features of cluster emis- 
sion and its relation to the phase transition behavior. 
The same "freeze-out" density of 0.38po is used to make 
a comparison for the results of I-CMD and I-LGM. Fig- 
ure 5 shows that N p , N n , N cp , Ni m f and A max changes 



with temperature within both models. The multiplicities 
of clusters and the largest fragment mass are close each 
other between I-LGM and I-CMD except for the mul- 
tiplicities of neutrons and protons, illustrating that the 
I-LGM is, in general, a good tool to describe the frag- 
mentation if Coulomb interaction can be ignored. When 
Coulomb interaction is switching on, N p , N cp and Ni m f 
increase due to the repulsive role among protons and 
while A max decreases. Meanwhile, N n does not change 
because of no Coulomb interaction. Ni m f also shows 
a rise and fall with temperature in the I-CMD cases. 
Coulomb force makes the turning temperature of Ni m f 
smaller due to the long range repulsion. 



The slopes of multiplicities of emitted clusters and of 
mean mass of the largest fragment have been plotted as a 
function of temperature in case of I-CMD in figure 6. The 
definite peaks of slopes are found as in the I-LGM case. 
The corresponding temperature at peaks locates about 4 
MeV in the I-CMD case with Coulomb interaction, and 
about 5 MeV in the I-CMD case without Coulomb inter- 
action. This turning temperature should also reflect the 
onset of phase transition there. If we investigate r, 52 
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FIG. 5: Same as figure 1, but for the comparison between differ- 
ent calculations: I-LGM (open circles), I-CMD without Coulomb 
interaction (open squares) and I-CMD with Coulomb force (solid 
squares). The "freeze-out" density of system is 0.38po- 



FIG. 6: Same as figure 2, but for the comparison between differ- 
ent calculations: I-LGM (open circles), I-CMD without Coulomb 
interaction (open squares) and I-CMD with Coulomb force (solid 
squares). The "freeze-out" density of system is 0.38po- 



and H: Figure 7, we will find that there are a minimum of 
r and the maxima of S2 and H around phase transition 
temperatures, i.e. about 4 - 4.25 Mev for I-CMD with 
Coulomb, around 5 MeV for I-CMD without Coulomb 
and around 5.5 MeV for I-LGM, respectively. Obviously, 
these phase transition temperatures are consistent with 
the ones extracted from the choppy position of the slopes 
of figure 6. 

Finally the RMS widths of the multiplicity distribu- 
tions of the clusters and of the largest fragment mass are 
checked in the I-CMD case in figure 8. The widths of the 
multiplicity distributions of neutrons and protons tend 
to be saturate at higher temperature, and while the ones 
for CP, IMF and A max demonstrate peaks at a certain 
fixed temperature, i.e. around 4 MeV for the case I- 
CMD with Coulomb, 5 MeV for the case of I-CMD with- 
out the Coulomb, which is similar to the I-LGM case. 
These turning temperatures are also consistent to the 
phase transition temperature as shown in figure 6 and 7 
in the I-CMD cases. Again, the largest fluctuation simul- 
taneously appears in the point of phase transition. 



Overall, the phase transition temperature seems to rely 
on some ingredients: such as the "freeze-out" density 
(or pressure), the model and its interaction potential, 
but the rule that emission rate and fluctuation of cluster 
multiplicity as a probe of phase transition is not changed. 



IV. CONCLUSIONS 

In conclusion, features of the emissions of light parti- 
cles, charged particles, intermediate mass fragments and 
the largest fragments are investigated for a medium size 
nucleus 129 Xe in the frameworks of isospin dependent lat- 
tice gas model and classical molecular dynamics model. 
N n , N p , N cp and A max show monotonously increasing or 
decreasing but Ni m f shows a rise and fall with temper- 
ature. Slopes of these observables go through extrema 
at the same temperature where the largest fluctuation of 
cluster multiplicity distributions is also observed. This 
temperature is consistent to the phase transition tem- 
perature extracted from the extreme values of r, S2 and 
H. It gives an indication that the cluster emission rate 
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FIG. 7: Same as figure 3, but for the comparison between differ- 
ent calculations: I-LGM (open circles), I-CMD without Coulomb 
interaction (open squares) and I-CMD with Coulomb force (solid 
squares). The "freeze-out" density of system is 0.38po- 



can be taken as a probe of the phase transition of nuclei 
and while the largest fluctuation is simultaneously ac- 
companied when the onset of phase transition occurs. In 
addition, the systematical comparison of I-LGM and I- 
CMD shows that LGM is a rather good tool to study the 
nuclear disassembly when the system is not large where 
the Coulomb interaction can be ignored. In light of this 
work, we think that the experimental study about the 
cluster emission is rather meaningful, especially to mea- 
sure the excitation function of the multiplicities, their 
slopes and variances, from which some signals for phase 
transition could be found. 
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